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Abstract 



We study a d = 2 discrete solid-on-solid model of complete wetting of 
a rough substrate with random self-affine boundary, having roughness 

^ . exponent (g. A suitable transfer matrix approach allows to discuss 

^ ! adsorption isotherms, as well as geometrical and thermal fluctuations 

Q I of the interface. For (g < 1/2 the same wetting exponent "0 = 1/3 as for 

P;~ I flat substrate is obtained for the dependence of the coverage, 6*, on the 

T^ij- ■ chemical potential, h {9 ^ h~'^ for /i — > 0). The expected existence of 

0\ . a zero temperature fixed point, leading to -0 = Cs/(2 — Cs) for Cs > 1/2, 

^y ! is verified numerically in spite of an unexpected, very slow convergence 

g I to asymptotics. 
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1. Introduction 

Adsorption of a nearly saturated vapour on a rougti substrate lias botli fundamental and 
technological relevance. In the case of complete wetting [1] , which is the main concern of the 
present paper, one considers the approach to infinity of the thickness of the adsorbed liquid 
film in presence of the vapour phase, as the chemical potential tends to the coexistence value. 
Indicating by / such thickness and by h the chemical potential difference, one expects / oc h~^ . 
Once assumed that roughness has well-defined scaling properties, globally characterized by 
a fractal dimension or a roughness exponent, the power law describing complete wetting 
should depend on these dimensions [2-9] . For example, a self-affine surface is characterized 
by a roughness exponent Cs, which gives the average transverse height fiuctuation, \AS\, 
sampled on a portion of the surface with projection of linear size L, as \AS\ oc L^^.* The 
full elucidation of the scaling behaviours describing wetting and/or the crossover phenomena 
possibly affecting their detectability, is a hard problem, whose solution should also suggest 
useful new ways of probing roughness in experimental samples. 

In recent years there have been many attempts to solve such issues, mostly on the basis 
of relatively simple phenomenological models and scaling considerations [2-7]. Scaling laws 
describing complete wetting have been proposed for both fractal [2,3] and self-affine sub- 
strates [4-7] . Attempts to apply such results to the analysis of recent experiments concerning 
the adsorption of N2 on deposited Ag [3,10] or cyclohexane vapor on Si substrates [11] have 
also been made. Many ambiguities however remain in this respect, because, even the recog- 
nition of a fractal or self-affine character of the samples seems to be sometimes a matter of 
controversy [3-5] [10]. An important step towards a satisfactory theoretical control of wetting 
on rough substrates should be made by studying in detail statistical models incorporating 
all essential features of the problem. Being free of the ambiguities still affecting experiments 
interpretation, such a study would provide crucial tests of the existing conjectures and of 
their ranges of validity. 

In the case of fractal substrates a model investigation of this kind has been undertaken 
recently by the present authors [8] . Due to the strong limitations imposed by the presence of 
overhangs for a fractal boundary, the statistical treatment had to be limited to a mean field 
approximation. However, neglecting interface fiuctuations is not expected to lead to substan- 
tial modifications of the asymptotic results in that case. Indeed, as we also confirm below on 
the basis of our results, wetting of an extremely rough surface, like a fractal, can be considered, 
in comparison to a self-affine one, as a T = phenomenon [12]. On the other hand, thermal 
effects are still adequately described in a mean field treatment of non asymptotic features of 
adsorption isotherms, such as layering and capillary condensation phenomena. However, the 
computational complexity inherent even in such mean field description, limits considerably 
the possibility of detecting asymptotic laws, like the expected power law dependence of the 
coverage on chemical potential [8]. 

Self-affine substrates can be modelled by profiles without overhangs. Due to this sim- 
plification, they are relatively more accessible to statistical treatment. At least for moderate 
self-affine roughness, thermal fiuctuations of the interface are expected to play a key role in 
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determining wetting exponents, as it is the case for fiat substrates when the bulk dimension 
is lower than 3 [12]. Thus, an adequate treatment of self-affine substrates must describe 
correctly interface fluctuation phenomena. In the present paper we accomplish this task 
by a transfer matrix study of an interface model of adsorption on self-afRne boundary in 2 
bulk dimensions. The interface is treated in a solid-on-solid (SOS) approximation and takes 
the shapes of a directed self-avoiding walk on square lattice. Thermal interface fluctuations 
are controlled by surface tension and lead to self-afiine behavior with a roughness exponent 
( = 1/2 [9]. The boundary of the substrate has the same kind of shapes, with quenched fluc- 
tuations. Generalization of a method introduced by Mandelbrot [13] enables us to produce 
random boundaries with preassigned roughness exponents, < Cs < 1- The liquid in the 
wetting film is further supposed to feel a van der Waals type of potential from the substrate, 
and its thickness / is controlled by the chemical potential. 

Existing scaling and analytical approaches to this problem [5,14] have strong ties with the 
work of Fisher and Lipowsky [12], who first discussed the role of the interfacial ( in determin- 
ing exponents like t/j in the flat case. The basic idea emerging from such scaling arguments, 
which will be recalled in the next section, is that for complete wetting the asymptotics should 
be determined by the maximum between C and Cs- When C > C^, the t/j exponent is the same 
as in the flat case, and thermal fluctuations control the asymptotics. When C < C^, control is 
assumed by the quenched geometrical fluctuations, and the role of ( is taken by (s in the for- 
mulas. Even if replica based perturbative renormalization group calculations have confirmed 
this simple scheme [14] , and the fact that only (s determines scaling in the strong fluctuation 
regime, with C^ > C, a direct verification did not exist and is one of the goals of our study. 

The application of transfer matrix techniques to a model like the one studied here encoun- 
ters some complications which we explain below, and involves relatively heavy computing. 
Thus, the approach presented in this paper has some methodological interest and should also 
constitute a starting point in view of further applications. 

This paper is organized as follows. In the next section we introduce the model, discuss 
some technicalities of the transfer matrix approach, and present a simple derivation of the 
scaling results which are the object of our numerical investigation. Section 3 is devoted to a 
discussion of test applications to the case of fiat substrate. In section 4 we present the main 
results for substrates with different degrees of roughness, with particular emphasis on the 
basic scaling laws involved. The last section contains general conclusions and outlines some 
open perspectives. 



2. SOS Model, Transfer Matrix Approach and Scaling Arguments 

In the SOS approximation the separation between condensed and vapour phases of the 
wetting fiuid is described by a sharp interface, which fiuctuates without forming overhangs. 
In 2 bulk dimensions on a square lattice this interface can be represented by a single valued 
integer function Z^ of the abscissa x, and fiuctuations are controlled by a Hamiltonian of the 
form: 

H = 5^[||Z, -Z,+i| + t/,(Z,) + /iZ,] (1) 



where x also runs over integers. The first term in the sum of eq.(l) represents the energy 
cost for the interface steps due to surface tension. Indeed, e corresponds to the strength 
of the attractive nearest-neighbor energy between fiuid particles in a lattice-gas model, in 
which particles can be placed at the centers of each elementary square not already occupied 
by substrate atoms. The second term in the sum is the energy due to long range van der 
Waals interaction with the substrate and h is the deviation of chemical potential from the 
coexistence value. 

The substrate boundary can be represented by a single valued function S^. Upon aver- 
aging over randomness, the following scaling is expected to hold 



\S,-S^>\r^\x-x'\'^% (2) 

Cs being the surface roughness exponent, with values in the interval (0, 1). 

The interface heights can be written as Z^ = Sx + Zx where Zx is the local thickness of 
the wetting film, i.e. the distance of the SOS interface from the substrate. 

Ux{Zx) in eq.(l) represents the potential energy, due to the substrate, of a column of fiuid 
with height Zx ai x. The long range nature of the potential and numerical checks convinced 
us that the following local approximation for this function 

Ux{Zx) = -U^^'^+c/z^-\ (3) 

is adequate to reproduce the scaling properties of adsorption in the high coverage regime. In 
all calculations reported here the choice a — 4 has been made. This corresponds to a 1/r^ 
van der Waals interparticle potential. 

We thus rewrite Hamiltonian (1) in the form: 



w = -Et'i°' + E 



7:1 Zx — Zx+i + ^Sx\ -\ — ^— Y + hzx 

.^ Zx 



(4) 



where i^Sx = Sx — Sx+i are the substrate surface steps and L is the horizontal linear size of 
the system. Here and in the following we assume 1 < z < oo. 

We have set up a transfer matrix technique [9,15] for the solution of the problem. In our 
case, because of roughness, and according to Hamiltonian (4), we can associate to any x two 
transfer matrices, T_(_(x) and T_(a;), with elements 

[T+(x)],,,, = exp [-(^1^ -z' + ASx\ + -^ + hz)/T\ (5a) 

[T_(x)],,,, = exp \-C-\z -z' - ASx-x\ + ^ + hz)lT\ {5b) 

Zi Z- - 

where T is the temperature measured in units of energy. The first term in eq.(4) does not 
enter in the definition of (5a-6) since it does not affect adsorption properties. One can see 
that [T_|_(x)]2^;z' and [T_(x)];2^2/ represent the Boltzman weighing factor for a step of the SOS 
interface from (x, z) to {x + 1, z') or {x — 1, z'), respectively. 

A film of adsorbed liquid can be represented by a very long SOS interface with,e.g., 
its ends pinned by the surface (see Fig.l). The average distance of the interface from the 



substrate calculated in its middle point will give a measure of the thickness of the wetting 
film. The evaluation, at a given temperature, of such thickness as a function of chemical 
potential h gives the adsorption isotherms. 

Let us consider a SOS interface of length L. This interface can ideally be separated 
into two parts at the middle point Xm'- for example, the left part has the left-end pinned at 
{xm — L/2j 1) and the right-end at {xm, z) (see Fig.l). We begin by treating the two parts 
as mutually independent. Let us consider first the left one. While its left-end is pinned, as 
we move to the right, the interface becomes more and more free to fiuctuate. The position of 
the interface at Xm can be studied in terms of an occupation profile r]^{xm,z) which is the 
probability for the interface to be at (x^,, z). It is easy to verify that 

ri+{xm, z) oc [T+{x^ - L/2) ■ T+(x^ - L/2 + 1) ■ ■ ■T+(x^ - l)]i,,. (6) 

This follows from iterating the formula: 

r7+(x + l,z)= iV++i Y. '?+(^' ^')[T+(x)],^, (7) 

where N^^^^ = 1/ "^^z ^,, rj{x,z')[T^{x)]z',z" is a normalization constant C^^V+i^^ ^) — 1) 
and r]^{xm — L/2, z) = 6^,1- In general, the shape of the occupation profile in Xm will depend 
not only on L but also on the shape of the substrate surface between Xm — L/2 and Xm- 
However, we expect that for large L the interface will loose memory of the starting point, 
and r]+{xrmz) will depend only on the shape of the substrate within a distance ^y, to be 
defined below. 

For the right part of the interface we can obviously repeat similar considerations. The 
occupation profile in this case will be obtained by iterations of the form: 

r7_(x-l,^) = iV-_,^r7-(x,/)[T-(x)],,,, (8) 

z> 

The occupation profiles at Xm we can get from eqs. (7) and (8) will in general differ from 
each other (see Fig.l). This is due to the roughness of the substrate which makes the surface 
asymmetric. 

In order to determine a unique, continuous interface on the basis of r\^ and ?7_, it is 
natural to impose that the left and right interfaces match at x^n- We thus define a profile: 

rj{Xra,z) = N^^r]+{Xm,z)7]-{Xra,z) (9) 

where N^^ = 1/ J2z V+i^m, z)r]- {xm, z) is the corresponding normalization factor. Of course, 
the same procedure discussed for Xm and leading to eq.(9) can be applied to any x, leading 
to the construction of a unique profile, as illustrated in Fig. 1. 

To obtain macroscopic quantities, two average operations have to be performed. We 
will distinguish between thermal and disorder averages, marked by brakets and overbars, 
respectively. The main quantity we are interested in is the coverage, i.e. the thickness of the 
wetting film, 9. 9 is given by 



9={z{x)) = Yzr]{x,z) = Yzr]{x,z) (10) 



On the other hand, also interface roughness has to be considered. In the case of a flat 
substrate, interface roughness is due to thermal fluctuations alone. Here, we can not exclude 
that the interface roughness can also be affected by the self-affine substrate. Information on 
the interface is given by the difference correlation function 

1- 



AC(x - x') = -{[z{x) - z{x')]^) = (z^x)) - {z{x)z{x')) (11) 

The last thermal average on the r.h.s of eq.(ll) can be defined following lines similar to those 
discussed above for the case of a quantity depending on one x alone. An idea of how this can 
be achieved is given in the following section, where an application to the flat case is discussed. 
The roughness of the interface can be measured by the perpendicular correlation length ^j^ 
which is extracted from AC in the \x — x'\ — » oo limit 



ii = AC{oo) = {z^x)) - {z{x)) (12) 



It is interesting to note that, by summing and subtracting the quantity {z{x))'^ to eq.(12 
^± can be expressed in the form 

-,R 

with 



ei = elT + eli? (13) 



^It = {z^{x)) - (z(x))2 = (z^x)) - (z(x))2 (14a) 

and 



.2 ^—T^ ^—r-2 



a,R = {4x))^ - {z{x)) (146) 

Thus the width of the interface is the sum of a thermal term (^±,t) , which is just the quenched 
average of the local mean square fluctuation around the mean position of the interface, and a 
geometrical one, i.e. the mean square fluctuation, induced by the surface geometrical disorder, 
of the average mean position (^_L,i?). The second term clearly vanishes when a flat surface is 
considered. 



The parallel correlation length ^y can also be calculated from {z{x)z{x')). Indeed one 
expects [9] 

AC(oo) - AC(x - x') ~ exp(-|x - x'\/i\\) (15) 

in the |x — x'| ^ oo limit. 

As we verify in the next section, for an interface bound near a flat substrate, 

u ~ i\ (16) 

with C = 1/2 in (i = 2. This interface wandering exponent in (i = 2 is explained by the fact 
that a free interface performs a random walk in the vertical direction [9]. 

Following Lipowsky and Fisher [12], we can separate into two terms the free energy of a 
wetting fllm of average depth /. The flrst one is U{1) = hi + c/l'^~^ + const., and represents 
the free energy contribution due to the last two terms of the second sum in eq.(4). The 
second term can be written on the basis of the standard continuum interface Hamiltonian, 
which replaces the flrst term in the above mentioned sum by a gradient squared of the local 



depth. For our interface, this gradient can be represented by ^±/^\\ ^ C± ■ Moreover, 

since for a bound interface it is natural to expect / ~ ^_l, the free energy density due to 
interface roughness, takes the form: 

/, ^ t;/2(C-i)/C. (17) 

In the case of a flat substrate, and in the absence of any kind of randomness, for d = 2, 
( = 1/2 has to be assumed in eq.(17), and minimization of U + fj gives 

I — >oo 

for 0" > 3, which leads to / ~ 6* oc /i~^/^, i.e. ip = 1/3. 

If the substrate is rough, with exponent Cs, the problem arises to establish which ( 
should actually enter is eq.(17). It is rather natural to expect that, as long as Cs < 1/2, 
the intrinsic thermal roughness of the interface determines //. Thus C = 1/2 has to be 
chosen in eq.(17) and the flat case t(j = 1/3 still holds. On the other hand, for C^ > 1/2, 
the geometrical substrate roughness induces a corrugation energy for the interface, which can 
still be estimated via eq.(17), this time with ( = (g > 1/2. 

This energy term becomes now asymptotically dominant in the minimum condition: 

^(/,/_^'/2(C.-i)/C.) = o (19) 

(J/t 

yielding if; ~ Cs/(2 — Cs), Cs > 1/2. In eq.(19) we again exploited the fact that / ~ ^±. 

These results were first derived by Kardar and Indekeu [5] on the basis of similar scaling 
arguments. An important remark to be made concerns the fact that the replacement of 
C = 1/2 by C = Cs > 1/2 in eq.(17), leads to an // which violates hyperscaling. Indeed, 
since // is a free energy per unit horizontal length, it should scale like ^,7 ^ C± ■ For 
( = 1/2, taking into account that / ~ ^_l, /i? ~ /~^ is indeed consistent with hyperscaling. 
For C = Cs > 1/2, the roughness free energy in eq.(19) does not scale as ^,7 . The hyperscaling 
violation is consistent with a zero temperature fixed point describing complete wetting. In 
the regime Cs > 1/2, interface fiuctuations are of a quenched geometrical, rather than thermal 
nature. A direct verification of this will be obtained in section 4, where it will be shown that 
while Ci_ ~ C±,T » U,R for Cs < 1/2, U ^ U,R » U,t for Cs > 1/2. 



3. Flat Substrate 

As a preliminary test of the transfer matrix approach, we have considered the case of a 
flat substrate. The absence of geometrical disorder determines here several simpliflcations. 
In particular, the transfer matrices T^{x) and T_(x) in eqs. (5a-6) are equal to each other. 
Thus, a unique ^-independent matrix, T, has to be considered. As a consequence, the 
occupation proflles t]^ and ?]_ also coincide, and rj^{x,z) = r]_{x,z) = rjoiz). 



One should note that the present case could be treated, equivalently, by computing 
directly the partition function of the system, i.e. Z = Tr[T^]. For large sizes, L, of the 
system, Z can be approximated by the quantity A^ where A is the largest eigenvalue of T. 
The relation with our approach is due to the fact that rjo is the eigenvector of T associated 
to A [15]. 

Thus, the above defined occupation profile rj is now given by 'rj{z) = Nr]Q{z). In spite of 
the absence of disorder, rjo cannot be calculated analytically. So, we evaluated it numerically 
by iterating eq.(7) until convergence. The coverage 9 is then obtained by eq.(lO), in which 
obviously disorder averaging does not apply, ^j^ follows from eq.(12) or (14) (^_l,_r = 0). The 
limited complexity of the numerical calculations makes also the evaluation of the interface 
parallel correlation length, ^y, feasible. In fact, for the calculation of {z{x)z{x')), which is 
now equal to {z{x)z{x')), we can use the expression 

(zixjzix }) = —=^ , ,.^1 rr, , ,. , (20) 

E.,.'^W[Tl^-^'l].,.'r7(^') 

which straightforwardly follows from extending considerations made in the previous section. 

To work in complete wetting conditions the parameters e, c and T have to be conveniently 
fixed. In particular we have chosen the values c/e = 2, for which the system is wetted at any 
temperature (i.e. T^ = 0). Our choice is meant to represent a typical situation for complete 
wetting and is not far from values appropriate for real adsorption experiments, like krypton 
on graphite [16], argon on solid xenon [17], or N2 on Ag [3]. 

To test the asymptotics of adsorption, which corresponds to consider very thick adsorbed 
films, the size of T should be made as large as possible. Of course, in a numerical calculation a 
reasonable matrix size Zmax x Zmax has to be considered, realizing a compromise between the 
necessity of avoiding finite size effects and that of reducing computation times. In our runs 
Zmax ranged from a minimum of 300 (low coverage) to a maximum of 800 (high coverage) , in 
this flat case. 

As discussed in the previous section, due to the dominance of interfacial thermal fluctu- 
ations, adsorption isotherms should follow the law 9 ~ h~^'^ for /i ^ 0. On the other hand, 
at very low 6*, the effects of the attractive potential, U{1), possibly dominate with respect to 
//, and a behaviour I c::^ 9 ^^ h~^'" (see eq.(18)) could hold before the thermal fluctuation 
effects are able to impose the asymptotic regime. This -0 = l/cr is strictly a T = exponent, 
of the FHH type (see Ref. [18]), and the low coverage window in which it can manifest itself 
will be the narrower, the higher the temperature. 

In Fig. 2a) we report three isotherms at temperatures T/e = 0.2 (open squares), 
T/e = 0.4 (open circles) and T/e = 0.8 (heavy circles). In the isotherm at higher T the 
interface fluctuations are dominating already at very low 9, and the log-log plot of 9 ver- 
sus h (heavy circles) is everywhere consistent with t(j = 1/3 (slope of the dashed line). To 
further confirmation of this scaling behaviour, a fit of the isotherm with a function of the 
type Ah~^ + B yielded V' = 0.327 ± 0.004. For this case, the plot in the insert shows also 
that ^_L follows the same scaling law and essentially ^± c:^ 9 which is consistent with the 
assumptions made at the end of the previous sections. In Fig. 26), for the two cases at higher 
temperatures reported in Fig. 2a), the corresponding ^± are plotted versus the respective ^y 



obtained by eq.(20). The log-log plot shows that eq.(16) is obeyed with C = 1/2 (dashed line 
slope), as expected. Fitting of the data for the higher temperature case, with a function of 
the type ^^y^ resulted in a z/ = 0.478 ±0.002. The two isotherms represented by open circles 
and squares in Fig. 2a) demonstrate that, the lower is the temperature, the wider is the 
preasymptotic region of 6 for which the scaling behaviour is not far from the law 9 ~ h~^'^ 
((J = 4 in the actual example). The above results show that, at least in the case of flat 
substrate, complete wetting is correctly described by our transfer matrix approach. 



4. Complete Wetting of Self— AfRne Surfaces 

As mentioned in the introduction, our self-afRne boundaries are obtained by a random 
generalization of an algorithm introduced by Mandelbrot [13]. This consists in applying a 
recursive transformation to an initially staircase shaped lattice walk. In practice, given an 
even number n, we consider a starting directed walk of 2n^ steps, obtained by alternating 
elementary forward and upward steps on a 2-d square lattice. Then, given a second even 
number m, with m < n, the walk is divided in n equal parts and in {n — m)/2 of them, chosen 
at random, all the vertical steps are reversed. The same procedure is then applied to each of 
the n parts of the walk obtained after the first stage, and all proceeds for k times. It is easy 
to check that the profile obtained in this way has a roughness exponent (s = Inm/ Inn. We 
checked numerically that the average height-height correlation of relatively small samples of 
boundary profiles obtained as above satisfied the scaling law (2), with the expected Cs up to 
a percent, or so. 

Since our self-affine boundaries are random, quenched averaging is needed in our evalua- 
tion of the wetting film properties. In practice we were careful to use in each actual calculation 
big enough L's, so as to guarantee a high degree of self-averaging in each individual realiza- 
tion of the random profile. Another technical difference, compared to the fiat case, concerns 
the size of the transfer matrices. Due to the higher coverage induced by surface roughness, a 
larger Zmax has to be considered. Of course, Zmax increases with increasing C^. The results 
we obtained required to consider matrices with Zmax between 1000 and 1500, with lengths L 
up to 6 ■ 10^. 

To test the effects of surface roughness on adsorption properties, we have considered 
three values of Cs, Cs = 1/4, 1/2 and 3/4. According to the scaling arguments of the previous 
section, the first and third values are representative of the low and high roughness regimes, 
respectively. Cs = 1/2 marks the border line between the two regimes. For all Cs, occupation 
profiles were obtained by eqs. (7-9). The coverage 9 and ^j^, ^±^t and ^_l,_r were also 
computed, and, for 6*, the asymptotic scaling exponent i/j was extrapolated. Computing time 
limitations did not allow a direct evaluation of ^|| , at variance with the fiat case. The choice 
c/e = 2 discussed in the previous section was made in all cases. 

Fig. 3 reports the results for Cs = 1/4 boundaries. Temperature was fixed at T/e = 0.8. 
The coverage 9 (heavy circles) follows in the whole range of h values (about 6 decades), 
a scaling behaviour consistent with t(j = 1/3 (see dashed line). As a confirmation, a fit 
encompassing the whole data range with a function of the form Ah~'^ + B gave if) — 0.319 ± 



0.008. Heavy squares represent data for ^j^, which clearly follows the same behaviour as 
6* (6* ~ ^x). In addition, the thermal part ^x^t (open squares) gives the entire ^±, the 
geometrical one being much lower (open rhombs). Thus, as expected, interface roughness is 
dominated by thermal fluctuations. 

Fig. 4 shows the results for the case of boundaries with (s = 1/2, again at T/e = 0.8. 
In spite of the higher roughness of the substrate, 9 (heavy circles) and ^_l (heavy squares) 
still follow the same scaling as for flat substrate. A flt of the type described above gave 
ifj = 0.327 ±0.002, still very close to the expected i/j = 1/3. In this case ^±^t and ^_l,_r follow 
essentially the same scaling law, but, still, ^±^r « C±,t- 

In the case of high substrate roughness, it turns out that the asymptotic scaling law 
expected on the basis of the T = flxed point mentioned in the second section, does not 
show up easily, even at prohibitively high coverages. For high coverages thermal fluctuations 
are still very important and seem to induce a behaviour close to the t/j = 1/3 law. In order to 
reduce the role of thermal fluctuations and to obtain evidence of a crossover to the expected 
ifj = Cs/(2 — Cs)? W6 had to reduce the temperature to T/e — 0.5 this time, and to consider 
L = 6 ■ 10'^ and Zmax = 1500. With such choices, somehow at the limit of our computational 
possibilities, one detects a sort of saturation of the coverage (heavy circles in Fig. 5) at very 
low h. This saturation is due to truncation of the transfer matrices. On the other hand, our 
numerical checks indicate that truncation effects are essentially absent up to /i > 6 ■ 10~^. 
An interesting feature of the isotherm is that, for h just above this limit, the slope has a 
clear increase, and, for about one decade, seems to be consistent with a, t/j = 0.632 ± 0.021. 
This estimate is reasonably close to the value t(j = Cs/(2 — Cs) = 0.6 implied by the scaling 
arguments. We interpret these results as an indication that, after a preasymptotic regime of 
lower effective '0, the system tends to follow the expected law at very high coverages. In the 
preasymptotic regime the slope of the log-log plot of 9 is definitely lower, and rather close to 
the flat case one (dashed line). 

Fig. 5 also shows that ^_l (heavy squares) follows the same law as 9. Now, however, 
interface roughness is almost completely due to geometry, rather than temperature, and the 
roles of ^j^^T and $,±^r are interchanged with respect to the previous cases. ^_l,t (open squares) 
is much lower than ^x,/^ (open rhombs). It also turns out that, in the preasymptotic regime, 
C±,T obeys rather closely the h~^'^ law, as in the flat substrate case (dashed line). As 9 
reaches its expected asymptotic slope, also ^±^r rises more rapidly, and, before saturation, 
assumes a slope comparable with that of 9, as well. 



5. Concluding Remarks 

In this paper we performed a systematic study of the effects of self-affine substrate 
roughness on complete wetting properties of a fluid. This study was based on a discrete, 
d = 2 statistical model genuinely incorporating the essential features of the problem, namely 
interface fluctuations, long range substrate potential, and geometrical disorder of the surface. 

Up to now, existing studies of complete wetting on self-affine substrates all relied on 
essentially continuum descriptions, so that the problem of testing their scaling [4,5] on discrete 



in 



models remained open, together with that of precisely establishing the ranges of validity of 
various regimes and of locating their crossovers. 

We applied to our model a transfer matrix technique which allowed to make a useful 
distinction between thermally and geometrically driven interface fluctuations, and to relate 
them to the coverage properties. 

With a rather high level of accuracy we could verify that power law behaviours expected 
on the basis of scaling arguments are verified already at relatively low coverages, and for dif- 
ferent temperature choices, in the case of moderately rough substrates (Cs < 1/2). The case 
(s = 3/4 was expected to constitute an example of complete wetting exclusively controlled by 
geometrical disorder through a T = fixed point. As a matter of fact, a satisfactory treatment 
of this relatively high substrate roughness revealed a real computational challenge within our 
model. We could verify that with (s = 3/4 a very wide preasymptotic scaling regime with •0 
slightly larger than 1/3 exists. Only by suitably reducing the effects of thermal fiuctuations, 
responsible for the preasymptotic behaviour, it becomes possible to detect crossover to the 
expected t/j = 0.6, with an effort reaching almost the limits of our computational capabili- 
ties. Such slow set up of the asymptotic regime indicates the existence of important scaling 
corrections, and was not a priori expected. This should constitute a warning with respect to 
attempts to interpret numerical or experimental results in the field. 

The proof that a T = fixed point indeed controls asymptotics in the case of high rough- 
ness is an indirect confirmation that studies of complete wetting on fractally rough substrates 
can be correct in predicting the asymptotics, even if disregarding thermal fluctuations of the 
interface. This applies to our previous study of adsorption on a fractal substrate by mean 
field methods [8]. 

The problem of complete wetting does not at all exhaust the range of applicability of 
methods like those we applied in this paper. The natural further step, already undertaken by 
us, is the study of the critical wetting transition [1,9,19], possibly occurring when in models 
of the kind discussed here the temperature approaches from below some critical value, T^. 
The role of self-affine roughness in this case can not be discussed in terms of simple scaling 
ideas, and more sophisticated, albeit approximated, renormalization group methods have been 
produced to elucidate it [14] . In this context the role of our model studies will reveal even more 
crucial, also because, we anticipate, surface roughness can lead to new features of wetting, at 
low temperature, which do not seem to be easily catched by continuum descriptions. 
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FIGURE CAPTIONS 

Fig. 1 Shape of a SOS interface on a rough self-afRne surface with (g = 1/2 (represented by 
the stepped curve bounding the shaded area). The heavy continuous hne represents 
the mean shape of the interface with the ends pinned at A and B on the substrate 
and obtained by applying at each x a matching condition of the form (9). Dashed and 
long-dashed lines represent the mean shapes of the right and left interfaces discussed 
in the text. In the upper plot, the occupation profiles rj, rj+, rj- at the middle point 
(see arrow) are shown. The calculation were performed with c/e = 2, T/e = 0.5 and 
h = 2.44 ■ 10-4. 

Fig. 2 Adsorption on a flat surface. Calculations were done with c/e = 2. a) Adsorption 
isotherms for T/e = 0.2 (open squares), T/e = 0.4 (open circles) and T/e = 0.8 
(heavy circles). Dashed and long dashed lines express a h~^'^ and h~^'^ scaling laws, 
respectively. In the insert ^±_ (heavy squares) for T/e = 0.8 is plotted versus h. Also 
in this case the dashed line expresses a h~^'^ scaling law. h) ^±_ versus ^|| for the two 

1 /2 

cases at higher temperature reported in a). Dashed line expresses a ^n scaling law. 

Fig. 3 Adsorption on a rough self-affine surface with C,s = 1/4. The calculation were done 
with c/e = 2 and T/e = 0.8. Heavy circles and squares represent the value of coverage 
9 and ^_l as functions of h, respectively. Open squares (partially hidden by heavy 
squares) and open rhombs represent ^_l,t and ^_l,_r values, respectively. The dashed 
line obeys the h~^/^ scaling law. 

Fig. 4 Adsorption on a rough self-afRne surface with C,s = 1/2. The choice c/e = 2 and 
T/e = 0.8 was made in this case. Heavy circles and squares represent 9 and ^_l as 
function of h, respectively. Open squares (partially hidden by heavy squares) and 
open rhombs represent ^_l,t and ^j,,^, respectively. The dashed line gives again the a 
/j-i/3 power law. 

Fig. 5 Adsorption on a rough self-afRne surface with (s = 3/4. The calculation were done 
with c/e = 0.5 and T/e = 0.5. Heavy circles and squares represent 9 and ^± as 
functions of /i, respectively. Open squares and open rhombs give ^^ t and C±,R values, 
respectively. The solid line results from fitting 9 with a function of the form Ah~'^ + B 
in the range 6 ■ 10"^ < /i < 3 ■ 10"^. This gives ^ = 0.632 ± 0.021. The dashed line 
gives the h~^'^ scaling law. 
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